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^ ■ Abstract 

Q^ ■ Nuclear effects in the spin-dependent structure function g\ of the deuteron are 

^| , studied in the kinematics of future experiments at CEBAF, [y < 3 GeV, Q 2 < 

2 GeV 2 ). The magnitude of nuclear effects is found to be significantly larger than 
the one occurring in deep inelastic scattering (y — > oo, Q 2 —* oo). A possibility to 
. measure the neutron structure functions in the CEBAF experiments with deuterium 

C \ is analysed. It is found that disregarding or improperly treating nuclear effects in 

the region of nucleon resonances would lead to the "extraction" of an unreliable 

i— h ■ function. A procedure aimed at correctly extracting the neutron structure function 

from the deuterium data is illustrated and conclusions about the experimental study 

of the Q 2 dependence of the Gerasimov-Drell-Hearn Sum Rule for the neutron are 

(^) ■ drawn. 
t-h ; 

I. Recently it has been proposed at CEBAF to experimentally study the spin- 
dependent structure function (SF) of the neutron g™, in a wide interval of energy v ( 
0.2 — 3 GeV) and momentum transfers Q 2 (0.15 — 2 GeV 2 ), using polarized deuterium and 
3 He targets. These experiments will shed light on a number of quantum chromodynamics 
(QCD) sum rules and will help to establish a connection between results predicted by low 
energy theorems (Q 2 — > 0) and perturbative QCD (Q 2 3> m 2 , m being the nucleon mass). 

Of particular interest is the Q 2 dependence of the Gerasimov-Drell-Hearn (GDH) Sum 
Rule for the neutron. 

The GDH Sum Rule, which has been derived in the real photon limit (Q 2 = 0) by 
^ ■ Gerasimov and Drell and Hearn ||, reads as follows: 

m 2 7 dv k 2 

~ (0I/2O) - ^3/2^)) = —T-, (1) 
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where k is the nucleon anomalous magnetic moment, v t h is the threshold energy of the 
pion photo-production, 01/2(3/2) is the absorption cross section for total helicity 1/2(3/2), 
a is the fine structure constant. The sum rule ([I]) can be generalized to the case of electron 
scattering by expressing the helicity cross sections, 01/2(3/2), through the spin- dependent 
SF gi(x,Q 2 ) and g 2 {x,Q 2 ), obtaining M 



2m 2 



IgDh(Q ) = ~7y2~ I dx 
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Since under usual experimental conditions the second term in eq. (|2[) is small, one expects 
that 

Igdh(Q 2 = 0) = ~ (3) 

i.e. a negative value for the first moment of the nucleon SF. Since the experimental data 
in the region of the deep inelastic scattering show a positive value for the first moment 
of proton SF || [7|, ||, this implies that in the (Revolution of the GDH Sum Rule from 
the real photon limit to the deep inelastic region, the first moment of the proton SF must 
change its sign. The change in sign is expected to occur in the region of excitation of the 
nucleon resonances where neither perturbative QCD nor chiral theories are applicable. 
Indeed, the momentum transfer is already small enough to use the perturbative methods 
of QCD, whereas, at the same time, the chiral loops expansion does not work because of 
drastic changes of the helicity structure in the resonance region. 

On the other hand side, the measurement of the neutron SF will essentially contribute 
to the analysis of the deep inelastic sum rules, such as the Bjorken sum rule (BSR) and 
Ellis- Jaffe sum rules (see, for instance ref. || and references therein). In the deep inelastic 
limit, Q 2 —>■ oo, the BSR connects the difference of the first moments of the spin-dependent 
SF of the proton and the neutron with the axial constant of the neutron /9-decay: 



J dx {gl{x) - g^x)) = U^Y (4) 



where x = Q 2 /2mv is the Bjorken scaling variable. The experimental check of eq. (Q) has 
already begun || [7|, || by a measurement of the SF of the proton and the neutron, using 
in the latter case nuclear targets, viz. deuterium and 3 He. The experimental data from 
different groups give slightly different values for the BSR. However, theoretical efforts 
in computing the Q 2 -corrections recoincile the data with the theoretical prediction ([|), 
leading to the conclusion that the BSR is experimentally confirmed with an accuracy 
of about two standard deviations in the measured interval of Q 2 ||. Nevertheless, the 
investigation of the (Revolution of the BSR in a wide interval of Q 2 and the problem 
of its explanation within QCD remain of great interests. Beside the GDH and BSR sum 
rules, yelding integral characteristics of the SF, a detailed study of the resonance behavior 
of the nucleon SF g±(x, Q 2 ) is planned as well. 

All above examples demonstrate that the measurement of the spin-dependent neutron 
SF at CEBAF will provide us with important new information about the nucleon structure 
and will help to test a number of theoretical models and methods. 

Keeping in mind the lessons we have learned from the EMC-effect, one might expect 
that nuclear corrections could play an important role in estimating the neutron SF from 
the combined nuclear and proton data 0, 10[. In the region of finite Q 2 ~ m 2 and v ~ m, 



nuclear corrections are much more important than in the deep inelastic limit [0| . In this 
paper the role of nuclear structure effects in electron-deuteron scattering in the resonance 
region will be discussed, paying special attention to the procedure of the extraction of the 
neutron SF from the deuteron data in the kinematics of future experiments at CEBAF. 



2 



II. The nucleoli contribution to the deuteron structure functions is usually calculated 
by weighting the amplitude of electron scattering on the nucleon with the wave function 



of nucleon in the deuteron (for recent developments see e.g. [|12|, |13, [14], [15| and references 
therein). For the spin- dependent SF the most important effects are the Fermi motion and 
the depolarizing effect of the D-wave. Additional effects, such as off-mass- shell effects or 
nucleon deformation, are found to be small |T(| 17]. For finite values of Q 2 and v, the 
deuteron SF g®(x,Q 2 ) reads as follows 

d 3 k mv N ( „ 2 \ / £(x,Q 2 )k 3 



(2tt) 3 kq 

Vmax (x,Q 2 ) 



-Si(x/y,Q 2 )f D (y,Ux,Q 2 )), (6) 

y 

where = (gf + g™)/2 is the isoscalar nucleon SF and tyjj(k) the deuteron wave function 
with spin projection M. In the rest-frame of the deuteron, with q opposite the z-axis, 
kinematical variables are defined as: 

kq = u(k + £(x, Q 2 )k 3 ), k = m + e D - k 2 /2m, (7) 
f = q 3 /u = \q\/u = y/l + 4m 2 x 2 /Q 2 } Q 2 = -q 2 , x* = Q 2 /2kq, (8) 

where = —2.2246 MeV is the deuteron binding energy. The limits y m in(max)(x,Q 2 ) 
are defined to provide an integration over the physical region of momentum in (^) and 
to take into account the pion production threshold in the virtual photon-virtual nucleon 
scattering^. Since both y m in(%, Q 2 ) and y max (x, Q 2 ) are solutions of a transcendent equa- 
tion, explicit expressions for them cannot be given. However, in our numeric calculations 
they are accurately taken into account. 

Eqs. ©-(§[) have the correct limit in the deep inelastic kinematics (Q 2 — > oo, v — > 
oo). In this case: £(x, Q 2 ) — > 1, y min — > x, y m ax M D /m, and the usual convolution 
formula for the deuteron SF [ Efy |T5 is recovered: 



M D /m 

g?{x,Q 2 )= f %?( x / y ,Q 2 )f D (y). (9) 

Equation (||) defines spin-dependent "effective distribution of the nucleons", fo, which 
describes the bulk of the nuclear effects in gf. The main features of the distribution 
function, /d(|/), are a sharp maximum at y — 1 + eo/2m w 0.999 and a normalization 
given by (1 — 3/2Pd) (Pp being the weight of the D-wave in the deuteron). As a result in 
the region of medium values of x ~ 0.2 — 0.6 the deuteron SF gf(x) is slightly suppressed 
by a depolarization factor, (1 — 3/2Pd) xg±(x), compared to the free nucleon SF. However, 

2 For x not too close to the limit of single- nucleon kinematics, x — > 1, the quasi elastic contribution 
can be disregarded 
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the magnitude of this suppression is small (~ 1%) and this is why it is phenomenologically 
acceptable to extract the neutron SF from the deuteron and proton data by making use 
of the following approximate formula: 

g?{x, Q 2 ) « (l - \p d ) (fi{x, Q 2 ) + £{x, Q 2 ))/2. (10) 

In addition, when integrated over x, eqs. (0) and ([TOP give exactly the same result 
(r = Jdxg^x)), i.e. 



T D (Q 2 ) = (l - -P D j (T n (Q 2 ) + T p (Q 2 ))/2, (11) 

which allows to define exactly the integral of the neutron SF T n from the deuteron and 
proton integrals, without solving ([]). 

Eqs. ©-© at finite values of Q 2 and v are more sophisticated than the corresponding 
equations in the deep inelastic limit. In particular, they do not represent a "convolution 
formula" in the usual sense, since the effective distribution function fr> and the integration 
limits are also functions of x. This circumstance immediately leads to the conclusion that, 
in principle, when integrals of the SF are considered, the effective distribution can not be 



integrated out to get the factor similar to (1 — 3/2Pd) in (11). Another interesting feature 
of formulae (§])-(§) is the Q 2 -dependence of fjj and y m in,(max)(x,Q 2 ). If we again limit 
ourselves to the discussion of the integrals of SF, one concludes that the Q 2 - dependence 
of such an integral is governed by both the QCD-evolution of the nucleon SF and the 
kinematical Q 2 -dependence of the effective distribution of nucleons. 

Thus, we have established that in the non-asymptotic regime, equation (|TT|), in prin- 



ciple, does not hold. Furthermore, it is not clear whether an equation similar to (10 
could be applied in this region. Indeed, we are discussing the kinematical conditions 
pertaining to nucleon resonances, where the "elementary" nucleon SF explicitly exhibits 
Breit-Wigner resonance structures corresponding to the excitations of the nucleon by the 
photon and one expects that the Fermi motion and binding of nucleons will result in a shift 
and a smearing of the resonance structures. However, one can hope that actual effects 
will be quantitatively small so that eq. similar to ([11]) and ( |T0D could phenomenologically 
still be valid. 

III. In our numerical estimates we use a reliable parametrization of the proton and 
neutron SF given by Burkert ||18|| , which takes into account several nucleon excitations and 



provides a reasonable description of the available nucleon data in the resonance region. 



Using the Bonn potential model for the deuteron wave function []i~9 |, we carry out a 
realistic calculation of the deuteron SF, g®(x, Q 2 ) in the region of nucleon resonances. 

In Fig. |1] the results of the calculation of the deuteron SF, g±(x, Q 2 ) at Q 2 = 0.1 GeV 
and 1.0 GeV, are compared with the input of the calculation, i.e. the isoscalar nucleon SF, 
gi(x, Q 2 ). It can be seen that the role of nuclear effects in the resonance region is much 
larger (up to ~ 50% in the maxima of the resonances), than in the deep inelastic regime 



(~ 7 — 9%, depending upon the models |T2|, [13], [P], [T5[] , with resulting ~ 6 — 7% from 
the depolarization factor (1 — 3/2Pd) and ~ 1 — 2% from the binding effects and Fermi 
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motion). Such a drastic effect is a consequence of the presence of the narrow resonance 
peaks in the nucleon SF. 

Fig. |2| shows the results of extraction of the neutron SF from the deuteron and proton 
data by using the approximate formula ([10]) which, we believe to give an upper limit 
of the possible errors in this extraction. To emphasize the role of nuclear effects in the 
region of finite Q 2 , the extracted neutron SF is compared with the original (input in the 



calculation) parametrization of the neutron SF. The use of the approximate formula fllTf ) 
appears to be in some regions completely unreliable. This can be easily understood as 
follows: the proton and neutron SF have similar behavior in the resonance region, in 
that the positions of the nucleon resonances are the same for both of them, whereas the 
resonances in the resulting deuteron SF are smeared and shifted, compared to the isoscalar 
SF. Therefore, subtraction the proton SF from the deuteron one, in the maximum of 
the former, can result in a minimum for the neutron SF, instead of a maximum. The 
conclusion of our analysis is that nuclear effects in the resonance region are very specific 
and the approximate formula ( |T0"D does not work even for crude extraction of the neutron 
SF. Obviously, another method of extracting the neutron SF should be used. 

In ref. PUJ a rigorous method of solving eq. ([5]) for the unknown neutron SF has been 
proposed and applied in the deep inelastic region. It has been shown that this method, 
which works for both spin-independent and spin-dependent SF, in principle allows one to 
extract the neutron SF exactly, requiring only the analyticity of SF. It can also be applied 
by a minor modification to the extraction of the SF at finite Q 2 , which is our present aim. 

The basic idea is to replace the integral equation ([|) by a set of linear algebraic 
equations, KQ N = Gd, where K is a square matrix (depending upon the deuteron model), 
Qn is a vector of the experimentally known deuteron SF and Qn is a vector of unknown 
solution. Changing the integration variable in (|6]), r = x/y, we get 

g?(x,Q 2 )= J dTg?(T,Q 2 )±f D (x/T,Z(x,Q 2 )), (12) 

where T m i n (x^Q ) x/y max (x^Q ), T max (x^Q ) x max {Q ^}/ymin{XiQ ) and x max (Q ) is 
defined by the pion production threshold in virtual photon-nucleon scattering. Let us 
assume that the deuteron SF has been measured experimentally in the interval (xi,^) 
and a reasonable parametrization for the SF is found in this interval. Then, dividing both 
intervals (xi,X2) and (T min ,T max ) into N small parts, one may write: 

9?(xi,Q 2 )~t29?(fj,Q 2 ) I -fo^/r^dr, i = l...N, (13) 

>=i i 3 T 

where fj = T m i n + h(j — 1/2) and h = {j max — r min )/N. Equation fll3|) is already explicitly 
of the form Qd = KQn, therefore the usual linear algebra methods can be applied to solve 
it. 

Note that the range of variation of r is larger than the one for x. Therefore, in 
principle, the SF of the deuteron experimentally known in the interval (£i,£2) contains 
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information about neutron SF in wider interval (for example, in deep inelastic regime 
Tmin ~ x/2 and r max = 1). However extracting information beyond the interval f min = x\ 



to f max = %2 is almost impossible in view of the structure of the kernel of eq. (|12D and the 



kinematical condition of planned experimental data We have to redefine the kernel 



of eq. Q12Q to incorporate new limits of integration f min = x\ and f max = x 2 [pO 

The procedure of solving eqs. (|5|) in the kinematical region of finite Q 2 and v will be 
presented elsewhere in details; here we only stress that the method works with a good 
accuracy. To check it, we calculated the deuteron SF by formula (^) with the nucleon 
SF gi(x, Q 2 ) from ref. |18|] and the deuteron wave function of the Bonn potential (19 



Then the obtained gf{x, Q 2 ) has been used as "experimental" data to calculate the vector 
Qd in (|13|); the matrix K has been calculated by using the same deuteron wave function. 
Equation ( |l3l) has been solved numerically for various "experimental" situations (changing 
the "measured" interval (xi,x 2 ), for different Q 2 , etc.). The obtained solution, i.e. the 
extracted neutron SF, has been compared point by point with the input to the calculation 
of g®. We found that method is stable and allows one to unfold the neutron SF with 
errors not larger than 10 -4 , which is much smaller than the expected experimental errors 
(note, that in the present paper we discuss only the spin dependent SF, nevertheless all 
results and conclusions are valid for unpolarized structure function F2(x,Q 2 ) as well.) 

IV. In this section we discuss the role of nuclear corrections in the analysis of the 
integrals of the SF, such as the GDH Sum Rule. 

A very important observation has been made in the deep inelastic limit, namely the 
exact formula ([]) and the approximate formula fllCf ) give the same result for the integral 
of the neutron structure function, g™(x,Q 2 3> m 2 ) (see eq. ([□]))• The applicability of 
the approximate formula in the deep inelastic region is based on the conservation of the 
norm of the distribution f(y) by the convolution formula @. This circumstance can not 
be immediately extended to the case of the resonanse region, since: (i) the covolution is 
broken in eq. (§) and (ii) the normalization of the function f(y, x, Q 2 ) is different from one 
of f{y). The integral of the distribution f(y,x,Q 2 ) represents the "effective number" of 
nucleons "seen" by the virtual photon in the process when the virtual photon is absorbed 
by the nucleon and at least one pion is produced in the final state (it is less than 1 at low 
Q 2 and x ► x max ). 

However, surprisingly enough, the use of the formula ( |TTD in the resonance region 
gives results numerically very close to the integration of the exact equation (§). This is 
a consequence of the smallness of the effects breaking the convolution in eq. (|5|). These 
effects can be accounted for by a new equation: 

r D (Q 2 ) = (l - \P D ) N eff (Q 2 )(T n (Q 2 ) + T p (Q 2 ))/2 (14) 

Eq. flTJ]) and the integral of eq. (§) represent the definition of the "effective number" 
N e ff(Q 2 ); the latter depends upon the form of the nucleon SF g±, and, since this is 
expected to strongly oscillates (see Fig. 1), even the sign of the correction can vary. For 
instance we obtain using the SF from |L8] , 



N eff (Q 2 = 0.1 GeV 2 ) = 1.02, N eff (Q 2 = 1.0 GeV 2 ) = 0.997, (15) 
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i.e. a rather small effect (+2% and —0.3% correspondly) . Therefore eq. (13) appears to 



be a reliable one for estimating the integrals of SF: setting n e ff(Q 2 ) = 1 does not lead to 
errors larger than 3% for Q 2 = 0.1 — 2.0 GeV 2 . 

V. In conclusion, we have shown that the effects of nuclear structure in the extraction 
of the neutron SF in the resonance region are much more important than in the deep 
inelastic scattering. We have explained how the correct neutron SF can be firmly extracted 
from the combined deuteron and proton data. At the same time, we have found that the 
integrals of the SF, such as the GDH Sum Rule, can be estimated with accuracy better 
than 3% by the simple formula fllTD which is also valid in deep inelastic region. 

VI. We would like to express our gratitude to V. Burkert, S.B. Gerasimov and O. 
Teryaev for enlightining discussions. One of us (L.K.) thanks the INFN Sezione di Perugia 
for hospitality during his visit, when most of this work was completed. 
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Figure captions: 



Figure 1. The spin dependent structure functions gi(x,Q 2 ) for two values of Q 2 . The 
deuteron SF (solid line) is compared with the isoscalar nucleon SF (dotted line) used as 
input in the calculation of eq. (|). 

Figure 2. The neutron SF (solid line) extracted by the approximate formula (eq. (0)) 
compared with the original parametrization (dashed line) used in the convolution for- 
mula (H|). 
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Figure 1: C. Ciofi degli Atti et al, The Neutron Spin Structure Function from. . . 
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Figure 2: C. Ciofi degli Atti et al, The Neutron Spin Structure Function from. . . 
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